
Calhoun 

imU'uiioiu! A Kim of die Nav.il Po<ic^rdduacc School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


1994-06 

Laser-Doppler velocimeter measurements in a 
cascade of controlled diffusion compressor 
blades at stall 

Ganaim Rickel, Humberto Javier 

Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/30858 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


http ://w w w. nps.edu/lEbrary 


Calhoun is the Naval Postgraduate Sc hoof's public access digital repository for 
research mate hate and institutional publications created by tire NPS community, 
Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed — and published — scholar^ author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 



NAVAL POSTGRADUATE 
SCHOOL 

Monterey, California 



THESIS 


LASER-DOPPLER VELOCIMETER MEASUREMENTS IN A 
CASCADE OF CONTROLLED DIFFUSION COMPRESSOR 
BLADES AT STALL 

by 

Humberto Javier Ganaim Rickel 
June, 1994 

Principal Advisor: Garth V. Hobson 


Approved for public release; distribution is unlimited. 


Thesis 

G1435012 





DUDLF' KWX L 
NAVA; : >" 

MONltRL. «A 


rr> ARY 

A 3CH00! 

iJwJ'tJ-iilOl 



REPORT DOCUMENTATION PAGE 


I. TITLE AND SUBTITLE 

LASER-DOPPLER VELOCIMETER MEASUREMENTS IN A CASCADE O! 
CONTROLLED DIFFUSION COMPRESSOR BLADES AT STALL 


r. PERFORMING ORGANIZATION NAMES(S) AND ADDRESS(ES) 


). SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 


the author and do not reflect the official policy or position of the I 


i. DISTRIBUTION/AVAILABILITY STATEMENT 


Unclassified 















Approved for public release; distribution is unlimited. 


LASER-DOPPLER VELOCIMETER MEASUREMENTS IN A CASCADE OF 
CONTROLLED DIFFUSION COMPRESSOR BLADES AT STALL 


by 

Humberto Javier Ganaim Rickel 
BS, Venezuelan Naval School, 1985 


Submitted in partial fulfillment 
of the requirements for the degree of 

MASTER OF SCIENCE IN ENGINEERING SCIENCE 

from the 


Author: 


Approved By: 


NAVAL POSTGRADUATE SCHOOL 



ABSTRACT 


An incipient compressor blade stall has been generated and examined in the Low 
Speed Cascade Wind Tunnel at the Turbopropulsion Laboratory. The test blades were a 
controlled-diffusion design with solidity 1.67, and stalling occured at 10 degrees of 
incidence above the design inlet air angle. Tufting and laser-sheet flow-visualization 
techniques showed that the stalling process was unsteady, and occurred over the whole 
cascade of 20 blades. Detailed laser-doppler velocimeter measurements over the suction 
side of the blades showed regions of continuous and intermittent reversed flow. The 
measurements of the continuous reversed flow region at the leading edge were the first 
data to be obtained of flow within the leading edge separation bubble. The intermittent 
reversed flow region measurements quantified what was observed in the flow 
visualization studies. Blade surface pressure measurements showed a decrease in normal 
force on the blade as would be expected at stall. 
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I. INTRODUCTION 


A. BACKGROUND 

The continuing effort to predict off-design performance and stalling behavior of 
compressor blades during the design phase has prompted studies to characterize the flow 
in and around leading edge separation bubbles of blades in cascade. Experimental studies 
have attempted to map viscous flow development in specific geometries. Recently 
Hobson and Shreeve [Ref. 1 ] reported detailed two-component (LDV) measurements of 
the flow through a controlled-diffusion (CD) compressor cascade at a Reynolds number 
of about 700,000, and at a very high-incidence angle (8 deg above design). 

They obtained a laminar leading-edge separation, which reattached turbulent, and 
enclosed a (laminar) bubble on the suction surface of the blade. Consistent with 
measurements at lower incidence angles, the reattached suction surface boundary layer 
remained turbulent and attached over the rear part of the blade. Since boundary layer 
separation ( for a code-validation test case) had not been achieved, the next step was to 
increase the incidence angle further to 10 deg above design and try to stall the (CD) 
blades. This was the motivation for the present study in which the flowfield through the 
CD cascade was extensively surveyed at a fixed incidence angle which was 2 deg greater 
than the previous incidence reported by Hobson and Shreeve [Ref. 1 ]. 


B. PURPOSE 

The objective was to increase the inlet air angle beyond 48 degrees, as tested by 
C'lassick [Ref. 2], Murray [Ref. 3], Hobson and Shreeve [Ref. 1], and Wakefield [Ref. 4], 
to 50 degrees in an attempt to stall the blades. The intention was to determine the 
maximum turning or lift generated by the blades, and to determine the way in which the 
suction-side boundary layer separated. Would the leading-edge separation bubble grow or 



would separation begin from the trailing edge where the boundary layer was fully 
turbulent. Two-dimensional laser-doppler velocimeter measurements were to be taken in 
the pitchwise or blade-to-blade direction at most of the stations measured by Hobson and 
Shreeve [Ref. 1], 

Prior to performing the above study, LDV measurements at 48 degrees were 
obtained in the inlet region in order to verify the results that both Hobson and Shreeve 
[Ref. 1] and Wakefield [Ref. 4] obtained during their experiments. This was desirable 
because Hobson and Shreeve had used different inlet guide vanes (lGV's) and, after new 
IGV's were installed, Wakefield performed only Hot-Wire measurements. A comparison 
of the measurements taken by the present author with those taken by Hobson and Shreeve 
at 48 degrees is presented in Appendix A. The study carried out at an inlet-air angle of 
50 degrees is reported in the sections which follow. 
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II TEST FACILITY and instrumentation 


A. LOW-SPEED CASCADE WIND TUNNEL 

The subsonic cascade wind tunnel and operating instrumentation were as described 
by Wakefield [Ref. 5], The cascade had 20 blades, the flow Reynolds number, based on 
chord length, was approximately 700,000 and the inlet air angle was 48 and 50 deg. 

The blades had a chord length of 5.01 in. and a spacing of 3 in. The blade 
coordinates and cascade geometry were reported by Elazar [Ref. 5], Figure 1 shows the 
schematic of the cascade. 

B. INSTRUMENTATION 

1. Pneumatic Data Acquisition System 

Blade surface static pressure measurements were recorded with a 48-channel 
Scanivalve. The pneumatic data acquisition system was the same as that described and 
used by Classick [Ref. 2] and the program "ACQUIRE" was used to perform the 
pressure measurements. Figure 2 shows the location of the pressure taps on blade number 
10, the location of which is shown in Figure 1. 

2. Laser-Doppler Velocimeter 

The horizontal (U) and vertical (V) mean velocity components, U-turbulence, 
V-turbulence, and Reynolds stress were measured with a two-dimensional LDV system 
consisting of four major subsystems:(a) the laser and optics, (b) the data acquisition 
system, (c) the automated traverse table, and (d) the seeding probe. A photograph of the 
LDV equipment, traverse table, counters and oscilloscope is shown in Figure 3, which 
also shows the north endwall of the cascade. 


























a. Laser and Optics 

A four beam, two color TSI model 9100-7 LDV system was used. The laser 
was a Lexell four-Watt Argon-Ion laser which was operated nominally at 2 Watts in a 
multi-line mode. Two colors, green (514.5 ran) and blue (488 nm) were selected by the 
color separator. The two beams were centered and split into a four beam arrangement to 
measure two velocity components at right angles to each other. Two Bragg cells shifted 
the frequency of one beam in each pair to allow measurement of reverse flows. The four 
beams then passed through a divergence section which improved the dimensions of the 
measuring volume. Two photo-detectors collected the scattered light after it passed 
through the same optics. Table 1 contains a summary of the characteristics of the LDV 
system. 


TABLE 1 


CHARACTERISTIC 

BLUE BEAM 

GREEN BEAM 

WAVELENGTH 

488 nanometers 

514.5 nanometers 

FRINGE SPACING 

4.51 microns 

4.76 microns 

FOCAL LENGTH 

762 mm 

762 mm 

NUMBER OF FRINGES 

28 

28 

HALF ANGLE 

3.10 degrees 

3.10 degrees 

MEASURING VOL. DIAM 

133 micro meter 

133 micro meter 

MEASURING VOL. LENG 

2.5 mm 

2.5 mm 

FREQ. SHIFT (FIND) 

+ 5 Mhz 

+ 5 Mhz 

BEAM SPACING 

82.5 mm 

82.5 mm 

ORIENTATION 

HORIZONTAL 

VERTICAL 

CHANNEL 

2 

1 

FREQUENCY SHIFT 

5 Mhz UP 

5 Mhz DOWN 


b. Data Acquisition 


The data acquisition system consisted of two TSI Model 1990 counter-type 
signal processors and a 1998A Master Interface in which the signals from the 
photo-detectors were digitized. An oscilloscope attached to the input conditioner of the 
counters provided real-time display of the photomultiplier output. The digitized signals 
from the counters were send to an IBM PC in which the information was processed by 





TSI proprietary software "FIND" version 4.0 . Through this software it was possible to 
position the LDV at programmed locations and automatically take measurements in 
surveys at any desired increment. 

c. Automated Traverse table 

The automated three-axis traverse was Model 9500 from TSI. The traverse 
used stepping motors for positioning the optical table which rested between the upper 
support arms. Digital encoders along each axis provided positioning accuracy to 0.0001 
inch. The traverse and encoder interface to the PC used RS-232C protocol. 

d. Atomizer and Seeding Probe 

Olive oil was used as a seed material in a TSI atomizer which produced 
approximately 1 micro-meter sized particles as measured by Elazar [Ref. 5], The seeding 
wand was adjustable, however, the adjustment was done on an arc, perpendicular to the 
tunnel, thus the seeding was not always at midspan. This limited the distance over which 
the pitchwise surveys could be extended. Figure 4 shows the atomizer and seeding probe. 






III. EXPERIMENTAL PROCEDURE 


A. PRESSURE MEASUREMENTS AND FLOW VISUALIZATION 

Once the tunnel was set up at 50 degrees and running at a plenum pressure of 12 
inches of water (approximately 700,000 Reynolds number), the pressure measurements 
were taken as specified by Classick [Ref. 2]. 

The flow visualization was carried out by projecting a laser sheet from the bottom 
left of the cascade to blade number 14, and while the tunnel plenum pressure was set at 
12 inches of water (gauge), fog was introduced through one of the endwalls. The flow 
pattern of the fog between the blades was illuminated by the laser sheet. This process was 
performed at night for better visibility. The process was filmed ussing an 8mm video 
camera. 

B. TUNNEL SET-UP AND TEST-SECTION CONFIGURATION 

For the present study, the 50 degree inlet flow angle was set by adjusting the inlet 
guide vanes and side walls to equalize the endwall static pressures on both upstream 
walls. The exit flow angle was adjusted by setting the tailboards at angles which gave 
nearly uniform downstream wall static pressure measurements in the pitchwise direction 
across the cascade. The average inlet flow angle was measured, with the LDV, over three 
passage widths, 31.3% of an axial chord length upstream of the blade leading edge. Fine 
adjustments of the inlet guide vanes were made to achieve an average inlet flow angle (as 
measured by the LDV) of 50.21 degrees. 

Previous LDV measurements were taken between blades 7 and 8 which were 
anodized black to minimize reflections. Because of the present inlet flow angle setting of 
50.21 deg., blade 8 was too close to the edge of the window. Thus blade 8 and 6 were 
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exchanged and all subsequent measurements during this study were taken between blades 
6 and 7 as shown in Figure 5. 

The tunnel reference velocity (Vref) was determined using the analysis of Elazar 
[Ref. 5]. At different tunnel speeds, the inlet flow velocity was measured (31,3% axial 
chord upstream) with the LDV, and the plenum pressure and temperature and ambient 
pressure were recorded. A least-squares curve fit was applied to the data to determine the 
calibration curve. During each subsequent run, the plenum and atmospheric conditions 
were recorded and used as input to a Newton method iteration algorithm to determine 
Vref. The result of this calibration is presented in Appendix E. 

C. LASER SET-UP 

The green beams of the laser were aligned vertically with the unshifted beam at the 
bottom and the blue beams were horizontal with the unshifted beam to the right , as 
shown in Figure 6. All surveys were conducted with the LDV optics "standard", i.e., the 
488-nm blue beam measuring the horizontal velocity component (U), and the 514.5-nm 
green beam measuring the vertical velocity component (V). Down shifting was used in 
the following form; the green beam was downshifted by 5MHz and the blue beam was 
upshifted by 5MFIz. The 1990 signal processors had the following settings: continuous 
(CONT) Mode; High Filter, 20MHz; Low Filter, 0.3MHz; Amplitude Limit, full 
counterclockwise; Cycles/Burst, 8; Comparison, 1 percent; Auto (green button), in; 
Voltage, External (EXT); Data Ready, lntemal(lNT); Gain, One (01); Resolution 
(No/SEC), One (01). For the Data Interface Master; Coincidence Mode, Range XI and 
Delta Interval 2 to the power 3 micro-seconds was used throughout this study. 

In the Optics screen of the acquisition menu of FIND the frequency shift was set to 
+5MHz on both channels. As the maximum reverse flow Doppler frequency was 
approximately 1MHz this level of 5MHz downshifting allowed the determination of 
reverse flow velocities, both in the mean and intermittently. The determination of this 
final selection is shown in Appendix D. 
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Figure 6. LDV Fringe Pattern and Beam Arrangement 












D. SURVEYS 


1. Inlet Surveys at 48 and 50 Degrees 

All LDV measurements presented herein were averaged over 3000 data points, 
and plus or minus 2 Standard Deviation histogram editing was performed for the 
flowfield distribution plots. The edited histograms were used to determine the edge of the 
separation and reverse flow regions. 

The initial pitchwise survey at station 1 (Figure 7) was conducted over three 
passage widths to determine the flow periodicity. All subsequent inlet pitchwise surveys 
were traversed over a 4 in. distance, spanning the region of maximum seeding. The first 
three inlet surveys, at stations 1, la and lb, were carried out with the LDV horizontal. 
Station lb was repeated with the laser pitched upwards by 4 deg. The need for pitching 
was to allow for closer access to the leading edge, i.e., so that there would not be any 
blade shadow interference at the subsequent stations lc-le. At any time during the 
experiment, if the laser was either pitched or yawed, then the previous survey would be 
repeated to enable the determination of any errors due to the measurement volume 
orientation. The maximum spatial error, due to probe volume orientation, was calculated 
by Hobson and Shreeve [Ref. 1] to be 0.3mm. This error was because the probe volume 
was not parallel to the blade span, and therefore seed particles displaced from the actual 
measurement location could be measured. The location of the measurement volume was 
always referenced to the same location between the blades throughout the study. The 
alignment procedure is described by Elazar [Ref. 5], 

2. Passage Surveys at 50 Degrees 

Measurements were taken only on the suction side, over a two inch pitchwise 
distance. Figure 7 shows the positions for the passage surveys and each dot on the figure 
represents a measurement location. These dots were stretched away from the surface to 
approximate a boundary layer survey. The passage surveys (between blades 6 and 7) were 
conducted with the same LDV optics configuration specified for the inlet surveys. In 
addition, the LDV was yawed by 4 deg to the left and pitched upward by 2 deg to avoid 
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Figure 8. Suction Side Passage Survey Locations 




the laser beams being shadowed by the blade. This was done for the suction side close to 
the leading edge, from station 2 to 7. At stations 7 to 15 the LDV was only yawed by 4 
deg. 

3. Wake Surveys at 50 Degrees 

Wake surveys (between blades 6 and 7) were conducted with the same LDV 
optics configuration specified for the inlet surveys. The LDV was horizontal and 
perpendicular to the tunnel for stations 16 to 19 and the surveys were performed over two 
passage widths (6 inches). Figure 7 shows the positions for the wake surveys. 
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IV. RESULT AND DISCUSSION 


A. BLADE SURFACE PRESSURE MEASUREMENTS 

The upper plot of Figure 9 shows the blade surface pressure distribution measured 
by Dreon [Ref. 6] at 40 and 43 degrees, Armstrong [Ref. 7] at 48 degrees and the present 
work at 50 degrees. The integration of the area within the pressure distributions for each 
angle gave the Normal Force Coefficient. The lower plot (Normal Force Coefficient 
versus Angle of Attack) shows a drop-off in force (or lift) at 50 degrees, consistent with 
the observation that the cascade had entered into stall. 

B. INLET SURVEYS (STATIONS 1 THROUGH IE) 

Figures 10 through 15 show the horizontal (U), vertical (V) components and the 
total velocity (Utot) distributions in the pitchwise direction ahead of the blades. At station 
1, a disturbance in the total velocity profile is evident which is periodic and three inches 
apart. This disturbance corresponds to the spacing of the blades and thus the presence of 
the blades is now felt 30% of an axial chord ahead of the leading edges. This magnitude 
of upstream disturbance, was not evident at lower inlet air angles. 

Station 1A (Fig. 11) shows measurement anomalies on the U component which are 
due to imperfections in the acrylic window. In subsequent figures (12 through 15) the 
total velocity (Utot) decreased as the flow approached the leading edge of blade number 6 
and then increased again as the flow rounded the leading edge of the blade. 

The final inlet profile (Fig. 15) shows a variation in total velocity of 40% (from 1.0 
to 0.6) across the leading edge. This variation is less than that previously measured at 48 
degrees inlet air angle (> 50% variations), and this too is an indication that stall had 
occurred. 
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Figure 10. Survey at Station 
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Figure 11. Survey at Station 1A 
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Figure 12. Survey at Station IB 
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Figure 13. Survey at Station 1C 
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Figure 14. Survey at Station ID 
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Figure 15. Survey at Station IE 



C. PASSAGE SURVEYS (STATIONS 2 THROUGH 15) 

At station 2 only forward moving particles were measured, and the mean velocities 
(both U and V components) were all positive (Fig. 16). The discontinuity in the V/Vref 
profile between points 11 and 12 was unexplained. At station 2A the magnitude of the 
first data point dropped off significantly (Fig. 17). Upon examination of the histograms 
for the vertical velocity component it was found to contain reverse flow particles, which 
indicated that this region had intermittent reverse flow. The first data point at station 2B 
had a negative mean V velocity and a positive mean U velocity (Fig. 18), and this 
indicated the beginning of the leading edge reverse flow region (i.e., negative mean 
velocity on V). The following 5 data points had intermittent reverse flow histograms. 

At station 3 the first three data points had negative mean velocities, both U and V, 
and then the following 7 data points had intermittent reverse flow particles. Station 4 only 
had intermittent reverse flow particles (no histograms with a negative mean) for the first 6 
data points. The discontinuity in the profile as shown in Figure 20 illustrated the change 
over from intermittent reverse flow to all positive, or forward-moving particles. The 
profile at station 5 (Fig. 21) was very similar to that at station 4. 

At station 6 (Fig. 22) the first two points showed only forward moving particles, the 
third data point had intermittent reverse flow, the next five data points were all positive 
and the ninth data point again had intermittent reverse flow. All other data points beyond 
the tenth point had histograms with only positive values. The first data point at station 7 
(Fig. 23) only had positive moving particles, the second through sixteenth data points 
showed intermittent reverse flow and then all higher points were positive. 

The first data point at station 8 (Fig. 24) had only positive particles, the next 17 
data points showed intermittent reverse flow, and then all the points showed only positive 
flow. The mean flow profile once again showed a significant discontinuity in that region. 

Stations 9 through 15 (Figs. 25 through 31) were similar in that they all showed 
regions of intermittent reverse flow close to the suction surface of the blade followed by 
the core flow where all the measured particles had positive velocity components. 
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Figure 16. Survey at Station 
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Figure 17. Survey at Station 2A 
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Figure 18. Survey at Station 2B 
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19. Survey at Station 
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Figure 20. Survey at Station 
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Figure 21. Survey at Station 
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Figure 22. Survey at Station 
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Figure 23. Survey at Station 
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Figure 24. Survey at Station 
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Figure 25. Survey at Station 
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Figure 26. Survey at Station 10 
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Figure 27. Survey at Station 11 
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Figure 28. Survey at Station 12 




Figure 29. Survey at Station 13 
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Figure 30. Survey at Station 14 
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Figure 31. Survey at Station 15 





It can be seen in Appendix B that for each station all the histograms to the left of 
the discontinuity had negative and positive velocities and the histograms to the right of 
the discontinuity had only positive velocities. 

D. WAKE SURVEYS (STATIONS 16 THROUGH 19) 

Figure 32 through 34 show the horizontal (U) and vertical (V) velocity components 
and the total velocity (Utot) distributions through the wakes at the exit of the cascade. 
Like the other surveys, each point in these plots represents a histogram of 3000 data 
points which where analyzed at plus and minus two standard deviations. The ones that 
delimited positive from negative velocities for each station are printed in Appendix C. 
Two features are evident in these plots; firstly, the width of the wake increased from 
station 16 to 19, and secondly, the region of intermittent reverse flow was within the 
wake on the suction side of the blade (to the left of the X=0 line for blade 6). 

E. SUMMARY 

Once all the histograms from each station were analyzed, the boundary of the 
region of intermittent reverse flow (last point of negative velocity at a station) was plotted 
for each station 2 through 19. This is shown in figure 36 with dotted lines. Also shown on 
this plot , with the solid line, is the region of reverse flow as determined by a negative 
mean velocity on the vertical component. This line represented the reverse flow region of 
the leading edge separation bubble which had been observed with flow visualization 
techniques. It was postulated that the reason reverse flow was measured in this region 
was because the flow was unsteady and seed particles were entrained into the leading 
edge separation bubble. This was not possible at lower inlet air angles because the flow 
was relatively steady compared to the present study. Flow visualization also confirmed 
the two distinct regions of intermittent reverse flow, as shown by the two regions of 
dotted lines; the lower region being associated with the leading edge separation bubble 
and the upper region representing the turbulent separation that occurred aft of mid chord. 
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Figure 32. Survey at Station 16 
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Figure 33. Survey at Station 17 
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Figure 34. Survey at Station 18 
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Figure 35. Survey at Station 19 
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Figure 36. Reverse Flow Regions 










More detailed surveys are needed between stations 6 and 7 to fully characterize the 
transition between these two regions. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The lack of experimental data of compressor cascades at or near stall has been 
somewhat alleviated with the current set of detailed measurements. The following 
specific conclusions can be drawn. 

1. The controlled diffusion (CD) cascade was successfully stalled. This was 
confirmed with the blade surface pressure measurements, which showed that for 
50 degrees the normal force on the blade had decreased. Flow visualization 
techniques (both tufting and laser sheet with fog or smoke) also confirmed that the 
blades had stalled. 

2. It was possible to measure both mean reverse flow and intermittent reverse flow 
with the LDV. With the appropriate use of frequency shifting it was possible to do 
these measurements with the certainty that the results from the histograms were 
correctly representing negative or positive velocities. 

3. The regions of reverse flow were plotted. With the information obtained from 
each histogram at each station it was possible to plot regions of intermittent 
reverse flow and also a region of leading-edge reverse flow. 

4. It was possible to take LDV measurements inside the reverse flow region during 
the stalling process, which was unsteady. 

5. The inlet-flow pitchwise surveys at an inlet air angle of 48 degrees compared very 
well with previous measurements. 
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B. RECOMMENDATIONS 

The following specific recommendations for further measurements at the -50 
degrees inlet-air angle setting are proposed; 

1. More detailed measurements should be taken in the leading edge separation 
bubble region (station 2 to 4). 

2. More detailed measurements should be taken between stations 6 and 8 to further 
characterize the region of forward and intermittent reverse flow. 

3. Detailed measurements are needed between stations 15 and 16 to determine the 
trailing edge base flow region. 

4. Pressure side passage surveys are also needed. 

5. Measurements away from mid span are needed to determine the degree of two 
dimensionality of the flow. 

Blade surface pressure measurements at approximately 49 degrees inlet air angle 
are also needed to determine the maximum blade loading condition. 
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Plfcchwise 




station 


Y(l) 


-4.896 
-4.896 
-4.896 
-4.896 
-4.896 
-4.896 
-4.896 
-4.896 
-4.896 


-4.896 

-4.896 

-4.896 

-4.896 


U/Vref 

1.032942 

0.899393 

0.611193 

0.592662 

0.596724 

0.609269 

0.626958 

0.65313 

0.685123 

0.726156 

0.78884 

0.884321 

1.019106 

0.846569 

0.627775 

0.609753 

0.616697 


V/Vref 

0.655777 

0.305473 

0.452963 

0.560462 

0.621573 

0.668145 

0.705764 

0.738446 

0.765836 

0.787 

0.799841 

0.784678 

0.655669 

0.34376 

0.48539 

0.585436 

0.645592 


IJ-turb 

refer 

1.95444 

4.000105 

1.218652 

1.266566 

1.264895 

1.407681 

1.471754 

1.490103 

1.481566 

1.542182 

1.883499 

1.922391 

2.272077 

1.513604 

1.529772 

1.495324 



2.395522 
4.453251 
1.492263 
1.320667 
1.247186 
1.275609 
1.314726 
1.267063 
1.287087 
1.306153 
1.530384 
2.078085 
3.191467 
2.621281 
1.485803 
1.327857 
1.342328 


Reynolds 
Stress 
0.309529 
-3.537 
0.128289 
0.138072 
0.166161 
0.15649 
0.242738 
0.210081 
0.121722 
0.117827 
-0.03132 
-0.77927 
-0.67551 
-0.66567 
0.192735 
0.26045 
0.165122 


Corre1. 
Coeff. 

0.09653 

-0.28991 

0.103002 

0.120522 

0.15379 

0.127247 

0.183168 

0.162463 

0.093201 

0.085408 

-0.01587 

-0.28482 

-0.13602 

-0.14038 

0.125133 

0.18721 

0.120114 
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B. HISTOGRAMS FROM STATION 2 THROUGH 15 FOR 50 DEG 
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Position (retain) (-1.1233, 0.0000, -4.7905) Velocity at cursor = 0.^12 

Velocity mean = 20.071 Percent at cursor = 0.00 
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Me locity (w/s) for Component 





















Velocity at cursor = 0/#* 
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Velocity (re's) for Component II 

Position (rec/ln) (0.7500, 0.0000, 0.6779) Velocity at cursor = - 0 .Z 08 

Velocity mean = 54.324 Percent at cursor = 0.00 



Position (recoin) (0.6250, 0.0000, 0.6779) Velocity at cursor = -0.008 

Velocity mean = 45.557 Percent at cursor = 0.23 
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D. TABLE OF SHIFT SELECTION AT PLUS OR MINUS 5MHZ AND LDV 


MEASUREMENTS. 


BLUE BEAM (NORMAL FLOW), FRINGES DIRECTION -► FLOW DIRECTION -» 



possible to measured a positive and negative velocity (normal and reverse flow). 
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E. TUNNEL CALIBRATION DATA 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccceccccoccc 

c OUTPUT FROM PROGRAM CALIBRATE C 

ccccccccccccccccccecccccccccceccccccccccccccccceccccccece(x.*cr(vcor<vccccrrcccr 

LEAST SQUARES STRAIGHT LINE CURVE FIT IS USED 
TO DETERMINE TUNNEL CHARACTERISTICS AT DIFFERENT SPEEDS 

NEWTON S METHOD IS USED TO DETERMINE THE REFERENCE VELOCITY 
FROM THE RECORDED AMBIENT PRESSURE AND TUNNEL PLENUM 
PRESSURE AND TEMPERATURE 


BEGIN DETERMINING TUNNEL CHARACTERISTICS 
FROM THE FOLLOWING MEASURED VALUES 


TANGENTIAL VEL. 
M PER SEC. 


AMBIENT PRESS. 
INCHES MERCURY 


19.9600 

32.0520 


24.6660 

39.2970 


29.8941 

29.8941 



PLENUM PRESS. 
INCHES WATER 


3000 

0000 

0000 


PLENUM TEMP. 


CALCULATED VALUES FOR THE TUNNEL CONFIGURATION 


TOTAL VALOCITY 


MACH NUMBER 


MACH NUMBER FUNCT. PRESSURE RATIO 


0.317303192E+02 
0.507107968E+02 
0.624438918E+02 
0.728529848E+02 
0.799598643E+02 
0.860418448E+02 


0.414931434E-01 
0.662819212E-01 
0.815012269E-01 
0.949967822E-01 
0.104164963E+00 
0.111981989E+00 


0.599998018E-02 
0.152081979E-01 
0.228644274E-01 
0.308775796E-01 
0.369544196E-01 
0.425268503E-01 


-0. 19331J650E * 03 
-0.83 685808 5E f-02 
-0.522360G85E^2 
-0.3797933OOE102 
-0.314827750Ef02 
-0.266449149F.I 02 


CALLING LEAST SQUARES SUBROUTINE 
TO DETERMINE THE PRESSURE RATIO AS A FUNCTION OF MACH NO. PARAM 

PRESSURE RATIO « A1 * ANUX + AO 


MATRIX EQUATION 

0.60E+01 -0.42E+03 AO 

-0.42E+03 0.50E+05 A1 

A1 » 0.19109469172E-03 


0.15E+00 
-0.71E+01 

AO = 0.39221597704E-01 


REFERENCE CONDITIONS FOR EACH RUN 
AMBIENT PRESSURE PLENUM PRESSURE 
INCHES MERCURY INCHES WATER 


PLENUM TEMPERATUF 
DEGREES CELSIUS 


RUM NAME 
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BEGIN NEWTON ITERATION 


ITERATION NUMBER 1 MACH NO. PARAM. 
ITERATION NUMBER 2 MACH NO. PARAM. 
ITERATION NUMBER 3 MACH NO. PARAM. 
ITERATION NUMBER 4 MACH NO. PARAM. 

VREF « 84.73722268314 


« 0.104363 ERROR TERM = -0.6401E-02 
- 0.110764 ERROR TERM - 0.1655E-03 
= 0.110598 ERROR TERM » 0.5999E-07 
= 0.110598 ERROR TERM = -0.1731E-10 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 



PARAM, 

PARAM. 

PARAM. 

PARAM. 


0.104165 ERROR TERM = 
0.110775 ERROR TERM = 
0,110598 ERROR TERM - 
0.110598 ERROR TERM - 
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PRESSURE RATIO = -31.35743 MACH NUMBER PARAMETER = 0.4151E-01 

RUN NAME = Q514slbl.PRN 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 


MACH NO. PARAM. 
MACH NO. PARAM. 
MACH NO. PARAM. 
MACH NO. PARAM. 


0.103675 

0.110804 

0.110598 

0.110598 


ERROR TERM * 
ERROR TERM « 
ERROR TERM = 
ERROR TERM = 




. 29953401308 


PRESSURE RATIO = -31.33544 MACH NUMBER PARAMETER = 0.4151E-01 

RUN NAME * 0514slcl.PRN 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 


1 MACH NO. PARAM. 

2 MACH NO. PARAM. 

3 MACH NO. PARAM. 

4 MACH NO. PARAM. 


- 0.103772 ERROR 

« 0.110798 ERROR 

* 0.110598 ERROR 

* 0.110598 ERROR 


TERM - -0.7025E-02 
TERM - 0.1994E-03 
TERM * 0.9887E-07 
TERM = -0.2852E-10 


VREF = 85.21942480550 


PRESSURE RATIO = -31.33544 MACH NUMBER PARAMETER = 0.4151E-01 

RUN NAME = 0514sldl.PRN 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 


MACH NO. PARAM. = 
MACH NO. PARAM. - 
MACH NO. PARAM. = 
MACH NO. PARAM. = 


0.103675 ERROR TERM « 
0.110804 ERROR TERM - 
0.110598 ERROR TERM = 
0.110598 ERROR TERM ■ 


85.29953401308 


-0.7129E-02 
0.2053E-03 
0.1066E-06 
-0.3074E-10 


PRESSURE RATIO = -31.33544 MACH NUMBER PARAMETER - 

RUN NAME * 0514slel.PRN 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 


1 MACH NO. PARAM. 

2 MACH NO. PARAM. 

3 MACH NO. PARAM. 

4 MACH NO. PARAM. 


= 0.103675 ERROR 
= 0.110804 ERROR 
** 0.110598 ERROR 
= 0.110598 ERROR 


TERM = -0.7129E-02 
TERM = 0.2053E-03 
TERM = 0.1066E-06 
TERM = -0.3074E-10 
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85.29953401: 


PRESSURE RATIO * -31.24322 MACH NUMBER PARAMETER = 0.4153E-01 

RUN NAME * 0511S21.PRN 


BEGIN NEWTON ITERATION 


ITERATION NUMBER X 
ITERATION NUMBER 2 
ITERATION NUMBER 3 
ITERATION NUMBER 4 


MACH NO. PARAM. 
MACH NO. PARAM. 
MACH NO. PARAM. 
MACH NO. PARAM. 


« 0.103968 ERROR 

- 0.110814 ERROR 

• 0.110625 ERROR 

* 0.110625 ERROR 


TERM * -0.6846E-02 
TERM = 0.189 3E-03 
TERM = 0.8621E-07 
TERM = -0.2489E-10 


VREF - 85.07919440271 


PRESSURE RATIO ~ -31.24322 MACH NUMBER PARAMETER - 0.4153E-01 

RUN NAME = 0511s2al.PRN 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 1 MACH 
ITERATION NUMBER 2 MACH 
ITERATION NUMBER 3 MACH 
ITERATION NUMBER 4 MACH 


PARAM. 

PARAM. 


0.103968 

0.110625 
0.110625 


ERROR TERM * 
ERROR TERM ** 
ERROR TERM = 
ERROR TERM * 


■0.6846E-C 
0.1893E-0 
0.8621E-C 


VREF = 85.07919440271 


PRESSURE RATIO = -31.24322 MACH NUMBER PARAMETER * 0.4153E-01 

RUN NAME - G511s2bl.PRN 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 1 
ITERATION NUMBER 2 
ITERATION NUMBER 3 
ITERATION NUMBER 4 


MACH NO. PARAM. 
MACH NO. PARAM. 
MACH NO. PARAM. 
MACH NO. PARAM. 


0.103772 ERROR 
0.110826 ERROR 
0.110625 ERROR 
0.110625 ERROR 


TERM = -0.7053E-02 
TERM - 0.2009E-03 
TERM = 0.1007E-06 
TERM - -0.2908E-10 


VREF = 85.23966280667 


I = 10 

PRESSURE RATIO * -30.93546 MACH NUMBER PARAMETER = 0.4155E-01 

RUN NAME = 0512S31.PRN 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 


MACH NO. PARAM. 


104066 


ERROR 


, 6770E-( 



ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 


MACH NO. PARAM. 
MACH NO. PARAM. 
MACH NO. PARAM. 


0.110836 

0.110651 

0.110651 


ERROR TERM 
ERROR TERM 
ERROR TERM 


VREF « 85.01903047416 


I ** 11 

PRESSURE RATIO » -30.95720 MACH NUMBER PARAMETER - 0.4155E-01 

RUN NAME - 0512S41.PRN 


BEGIN NEWTON ITERATION 

ITERATION NUMBER 1 MACH NO. PARAM. * 

ITERATION NUMBER 2 MACH NO. PARAM. = 

ITERATION NUMBER 3 MACH NO. PARAM. = 

ITERATION NUMBER 4 MACH NO. PARAM. = 

VREF - 85.09939011711 


0.103968 ERROR TERM = 
0.110842 ERROR TERM = 
0.110651 ERROR TERM - 
0.110651 ERROR TERM = 


-0.6874E-02 
0.1908E-03 
0.8790E-07 
-Q.2540E-10 


-31.22131 MACH NUMBER PARAMETER = 0.4153E-01 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 


MACH NO. PARAM. 
MACH NO. PARAM. 
MACH NO. PARAM. 
MACH NO. PARAM. 


0.103968 

0.110814 

0.110625 

0.110625 


ERROR TERM 
ERROR TERM 
ERROR TERM 
ERROR TERM 


‘0.6846E-02 
0.8621E“07 


15.07919440271 


PRESSURE RATIO = -31.22131 


MACH NUMBER PARAMETER 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 


1 MACH NO. PARAM. 

2 MACH NO. PARAM. 

3 MACH NO. PARAM. 

4 MACH NO. PARAM. 


= 0.103968 ERROR 

= 0.110814 ERROR 

* 0.110625 ERROR 

* 0.110625 ERROR 


TERM * -0.6846E-02 
TERM = 0.1893E-03 
TERM = 0.8621E-07 
TERM » -0.2489E-10 


VREF « 85.07919440271 


I = 14 

PRESSURE RATIO = -31.13385 MACH NUMBER PARAMETER = 0.4153E-01 

RUN NAME = 0507S71.PRN 
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ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 


0.103675 ERROR TERM - 
0.110859 ERROR TERM = 
0.110651 ERROR TERM = 
0.110651 ERROR TERM = 
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PRESSURE RATIO = -30.89209 MACH NUMBER PARAMETER = O.4155E-0] 

RUN NAME = 0507S111.PRN 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 


1 MACH NO. PARAM. 

2 MACH NO. PARAM. 

3 MACH NO. PARAM. 

4 MACH NO. PARAM. 


0.103578 

0.110865 

0.110651 

0.110651 


ERROR TERM 
ERROR TERM 
ERROR TERM 
ERROR TERM 


-0.7288E-02 
0.2145E-03 
0.1188E-06 
-0.3432E-10 


85.42010151240 


I = 19 

PRESSURE RATIO = -30.91383 MACH NUMBER PARAMETER = 0.4155E-01 

RUN NAME = 0507S121.PRN 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 


1 MACH NO. PARAM. 

2 MACH NO. PARAM. 

3 MACH NO. PARAM. 

4 MACH NO. PARAM. 


= 0.103578 ERROR 
* 0.110865 ERROR 
= 0.110651 ERROR 
= 0.110651 ERROR 


TERM « —0.7288E—02 
TERM = 0.2145E-03 
TERM = 0.1188E-06 
TERM = -0.3432E-10 


85.42010151240 


I = 20 

PRESSURE RATIO - -30.91383 MACH NUMBER PARAMETER = 0.4155E-01 

RUN NAME = 0507S131.PRN 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 


MACH NO. PARAM. 
MACH NO. PARAM. 
MACH NO. PARAM. 
MACH NO. PARAM. 


0.103578 

0.110865 

0.110651 

0.110651 


ERROR TERM 
ERROR TERM 
ERROR TERM 
ERROR TERM 




0.2145E-03 
0.1188E-06 
-0.3432E-10 


VREF » 85.42010151240 


I = 21 

PRESSURE RATIO = -31.51192 MACH NUMBER PARAMETER = 0.4151E-01 

RUN NAME - 0508S141.PRN 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 


MACH NO. PARAM. 
MACH NO. PARAM. 
MACH NO. PARAM. 
MACH NO. PARAM. 


0.103968 ERROR 
0.110786 ERROR 
0.110598 ERROR 
0.110598 ERROR 


TERM = -0.6818E-02 
TERM = 0.1878E-03 
TERM = 0.8454E-07 
TERM - ~0.2439E-10 
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VREF 


85.05899450070 


I = 22 

PRESSURE RATIO = -31.48982 MACH NUMBER PARAMETER « 0.4151E-01 

RUN NAME = 0508S151.PRN 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 


1 MACH NO. PARAM. 

2 MACH NO. PARAM. 

3 MACH NO. PARAM. 

4 MACH NO. PARAM. 


* 0.103870 ERROR 

« 0.110792 ERROR 

* 0.110598 ERROR 

* 0.110598 ERROR 


TERM - -0.6922E-02 
TERM - 0.1936E-03 
TERM = 0.9153E-07 
TERM « -0.2640E-10 


85.13925466046 


PRESSURE RATIO - -31.35743 MACH NUMBER PARAMETER = 0.4151E-01 

RUN NAME * 0514S161.PRN 


BEGIN NEWTON ITERATION 

ITERATION NUMBER 1 MACH NO. PARAM. 
ITERATION NUMBER 2 MACH NO. PARAM. 
ITERATION NUMBER 3 MACH NO. PARAM. 
ITERATION NUMBER 4 MACH NO. PARAM. 

VREF = 85.13925466046 


0.103870 
0.110792 
0.110598 
0.110598 


ERROR TERM = 
ERROR TERM = 
ERROR TERM - 
ERROR TERM = 


0.6922E-C 
0.1936E-C 
0.9153E-C 
0.2640E-1 


PRESSURE RATIO = -31.35743 MACH NUMBER PARAMETER - 0.4151E-01 

RUN NAME = 0514S171.PRN 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 


1 MACH NO. PARAM. 

2 MACH NO. PARAM. 

3 MACH NO. PARAM. 

4 MACH NO. PARAM. 


* 0.103772 ERROR 

- 0.110798 ERROR 

- 0.110598 ERROR 

= 0.110598 ERROR 


TERM * -0.7025E-02 
TERM * 0.1994E-Q3 
TERM « 0.9887E-07 
TERM - -0.2852E-10 


VREF = 85.21942480550 


I = 25 

PRESSURE RATIO = -31.35743 MACH NUMBER PARAMETER = 0.4151E-01 

RUN NAME « 0514S181.PRN 


BEGIN NEWTON ITERATION 

ITERATION NUMBER 1 MACH NO. PARAM. = 0.103772 ERROR TERM = -0.7025E-02 
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ITERATION NUMBER 2 MACH NO. PARAM. = 
ITERATION NUMBER 3 MACH NO. PARAM. - 
ITERATION NUMBER 4 MACH NO. PARAM. « 

VREF = 85.21942480550 


0.110798 ERROR TERM = 0.1994E-03 
0.110598 ERROR TERM - 0.9887E-07 
0.110598 ERROR TERM = -0.2852E-10 


BEGIN NEWTON ITERATION 


ITERATION NUMBER 1 MACH NO. PARAM. 
ITERATI0N NUMBER 2 MACH NO. PARAM. 
ITERATION NUMBER 3 MACH NO. PARAM. 
ITERATION NUMBER 4 MACH NO. PARAM. 


* 0.103772 ERROR TERM = -0.7025E-02 

- 0.110798 ERROR TERM = 0.1994E-03 

* 0.110598 ERROR TERM = 0.9887E-07 

« 0.110598 ERROR TERM = -0.2852E-10 











F. SURVEYS FROM STATION 1 THROUGH 19 
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